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Abstract 

We treat the classical dynamics of the hydrogen atom in perpendicular electric 
and magnetic fields as a celestial mechanics problem. By expressing the Hamiltonian 
in appropriate action-angle variables, we separate the different time scales of the 
motion. The method of averaging then allows us to reduce the system to two degrees 
of freedom, and to classify the most important periodic orbits. 

1 Introduction 

Our contribution to this Special Issue lies at the intersection of Martin Gutzwiller's sci- 
entific interests, namely celestial mechanics, electron motion, and chaos, especially its 
manifestations in quantal systems. We will be performing "celestial mechanics on a mi- 
croscopic scale" [ U&:91| by treating the dynamics of highly excited ( "Rydberg" ) electrons 
[C98] in crossed electric and magnetic helds using classical mechanics. 

Rydberg atoms in strong external fields constitute fundamental physical systems where 
the quantum mechanical regime of strong nonlinearity can be tested [|G90|, KvL9Sj. While 



the problem of a Rydberg atom interacting with a strong magnetic field (the Diamagnetic 
Kepler Problem, DKP, also known as the Quadratic Zeeman Effect, QZE) has been fairly 
well understood as a result of sustained research in the past two decades [ HRW89 , FrWi89| l , 
the superficially similar scenario resulting from the addition of a perpendicular electric field 

- the so-called crossed field arrangement fj|j| gSgg , |Wfc89| , |RFW91| , |vMDU94| , |vMU97a| 

- remains the least understood of all Rydberg problems. This is all the more remarkable 
in view of the prominence of the crossed fields in diverse areas of physics ranging from 
excitonic systems to plasmas and neutron stars. This problem is so complex because 
no continuous symmetry survives the extensive symmetry breaking DG89| ] induced by 
the two fields. The result is a wealth of new physics which is only possible beyond two 
degrees of freedom, such as Arnol'd diffusion [ TLL79j , p90| , LL92, vMDU96|. This absence 
of symmetry also allows localizing electronic wavepackets in all spatial dimensions, and 



the observation of these localized wavepackets [Y&93| has led to new insights into the 



'Dedicated to Martin C. Gutzwiller on the occasion of his 75th birthday. 

' Present address: Weierstrafi Institut, Mohrenstrafie 39, D- 101 17 Berlin, Germany 
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dynamics of the electron in the correspondence principle regime. It has also been found 
that a velocity-dependent, Coriolis-like force in Newton's equations causes the ionization of 
the electron to exhibit chaotic scattering [ MW92| , UF95 . JFU99]. All these phenomena, as 
well as renewed interest in the motional Stark effect [ JHY83| , F94|, make the crossed- fields 
problem an experimental accessible paradigm for a wide variety of outstanding issues in 
atomic and molecular physics, solid-state physics [DS92, 3chm93], nuclear physics JBM75 ], 
astrophysics | JMa89| |, and celestial mechanics |Mi82 ] 



The experimental challenge has been taken up by Raithel, Fauth, and Walther [ RFW91 



|RFW93| ] who in a landmark series of experiments have identified a class of quasi-Landau 
(QL) resonances in the spectra of rubidium Rydberg atoms in crossed electric and magnetic 
fields. Similar to the original QL resonances observed by Garton and Tomkins [GT69|, this 
set of resonances is associated with a rather small set of planar orbits of the crossed-fields 
Hamiltonian which is known to support an enormous number of mostly non-planar peri- 
odic motions |RFW91 ] . The dominance of planar orbits in these experiments has recently 



been explained [vMFU97a] 



In contrast with the DKP, and despite some preliminary work [FlWe9£], the systematics 
of periodic orbits in the crossed-fields problem has not been discovered up to now. The aim 
of the present work is to initiate a systematic classification of the orbits of the crossed-fields 
Hamiltonian, based on methods developed in celestial mechanics, specifically Delaunay 
variables and averaging. The analogy between atomic and planetary systems was already 



used by Born |B27|| , who studied in particular the crossed-fields problem, but neglected 
the quadratic Zeeman term because he was not studying Rydberg atoms, where it is 
prominent. 

Delaunay variables are action-angle variables which have a clear geometric interpre- 
tation in terms of Kepler ellipses. A fascinating historical account of the developments 
in celestial mechanics leading to the introduction of Delaunay variables can be found in 
a recent review by Gutzwiller [ |G98[ ]. These variables allow to separate the time scales of 
the motion, which is represented as a fast rotation of the planet (or the electron), along a 
Kepler ellipse with slowly changing orientation and eccentricity. 



The technique of averaging [V9€] allows to decrease the number of degrees of freedom 
by eliminating the fast motion along the Kepler ellipse. It has been used for a long time in 
celestial mechanics to compute the so-called secular motion of the solar system, and can be 
considered as a first order perturbation theory [ D69| , H70|. A systematic use of averaging 
allowed Laskar to integrate the motion of the solar system over several hundred million 
years fL9|, [L96|, |LR9|]. The method has been applied to the DKP in [pKN83| , |CDMW87|1 . 

We note in passing that the direct connection to celestial mechanics that the Delaunay 
variables provide is lost with an alternative set of variables called the Lissajous elements. 
These are obtained by regularizing the Coulomb Hamiltonian [F&92], and are appropriate 
for investigating the level structure of Rydberg atoms. The connection between the two 
sets in two dimensions is given by |DW91]. 

In this paper, we use the following notations. The Hamiltonian of an electron subjected 
to a Coulomb potential, a magnetic field B and an electric field F can be written in 
dimensionless units as 



H(Px,P y ,Pz]x,y,z) = -p 2 



l + | L . B + |(rAB)' 



r F, 



(1) 



where r = (x, y, z) is the electron's position, p = {p x ,Py,Pz 
its angular momentum. We write r = ||r|| and p = ||p||. 



its momentum, and L = r A p 
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Our paper is organized as follows. In Section ^, we summarize previous results on 
the case when only the magnetic field is present. This allows us to introduce Delaunay 
variables and the method of averaging, and illustrate them in a relatively simple situation. 

In Section H, we consider the case when only an electric field is present, which is 



integrable [R63J. We introduce another set of action-angle variables ("electric action- 



angle variables") based on parabolic variables, which are better adapted to perturbation 



theory [B27]. We then derive (new) transformation formulas from electric action-angle 
variables to the geometrically more transparent Delaunay variables. 

With these tools and sets of coordinates in mind, we finally turn to the crossed-fields 
Hamiltonian (jT]) in Section f|. We start by considering the two limiting cases B <C F and 
F <C B, which involve three distinct time scales, and can thus be analysed by a second 
averaging. We then study the dynamics in the plane perpendicular to B for general values 
of the fields. We conclude by an overview of the general structure of the phase space of 
the averaged Hamiltonian. 

2 The Quadratic Zeeman Effect (or DKP) 

We start by considering Hamiltonian (||) in the case F = 0. If we take the z-axis along 
the magnetic field B, (|l|) can be written as 

H = \ P *-\ + ^L z + ^ + yZ). (2) 

Although the equations of motion are easily written down, it is difficult to understand 
the qualitative properties of dynamics in cartesian coordinates. We will therefore take 
advantage of the fact that for small B, (||) is a small perturbation of the integrable Kepler 
problem, for which action-angle variables are known explicitly. By writing the Hamilto- 
nian (|2|) in these variables, we can separate slow and fast components of the motion. The 
qualitative dynamics can then be further analysed by using the method of averaging. 

2.1 Delaunay variables 

We start by considering the Kepler Hamiltonian 

1 9 1 

H = -p 2 --. 3 

2 r 

Besides the energy H, it admits as constants of motion the angular momentum L and the 
Runge-Lenz vector 

A = - + L A p. (4) 

r 

If H < 0, the motion takes place in a plane perpendicular to L, on an ellipse of eccentricity 
e = ||A|| and major axis parallel to A and of length 2a = —1/H. 

Action-angle variables taking these properties into account are well known in celestial 
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Figure 1. Definition of Dclaunay variables. The angles fi, i and u> determine the position 
of the Kepler ellipse in space. The line of nodes ON is the intersection of the plane of the 
ellipse with the xy-plane. The shaded area is proportional to the mean anomaly M . The 
true anomaly v and eccentric anomaly E are auxiliary quantities, allowing to relate the 
(X, y)-coordinates of the planet P with M. 

mechanics, where they are called Delaunay variables. The action variables are given byQ 

A = \fa 



G = vWl-e 2 = ||L|| (5) 



K = ^]~a\J 1 — e 2 cosi = L Zl 

where the inclination i G [0, tt] is the angle between L and the z-axis. These variables are 
denned on the domain \K\ ^ G ^ A. 

The corresponding angle variables are defined in the following way. The intersection 
between the plane of the orbit and the xy-plane is called line of nodes. The angle f2 
between line of nodes and x-axis is called longitude of nodes and is conjugated to K; the 
angle u between major axis of the Kepler ellipse and line of nodes is called argument of 
perihelion and is conjugated to G; the mean anomaly M, which is conjugated to A, is 
proportional to the area swept on the ellipse, according to Kepler's second law. 

In order to compute M, we introduce orthogonal coordinates (X, Y), where X is 
attached to the major axis of the ellipse. The true anomaly v and the eccentric anomaly 
E are auxiliary quantities defined by the relations 

X = r cos v = a(cos E — e) 

, (6) 

Y = r sin v = ay 1 — e 2 sin E, 

see Fig. ||. By trigonometry one can show that M and E are related by Kepler's equation 

M = E - es'mE, (7) 

which implies in particular that 

dM ^ r 

—— = 1 - ecos-E = -. 8) 
dE a w 



x We use the letters G and K in order to distinguish the action variables from their physical meaning. 
In celestial mechanics, one usually denotes L z by H instead of K, but we prefer the latter notation in 
order to avoid confusion with the Hamiltonian. 
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The transition from Delaunay variables to cartesian coordinates is done in the following 
way. Given the actions (|5|), we can compute 

G 2 K 
a = A 2 , e 2 = \-—, cos i = — . (9) 

We can then determine X, Y and the conjugated momenta 

1 sinE 1 



Px 



yfa 1 — e cos E 



(10) 







1 vi — e 2 cos E 1 
y 7 ^ 1 — e cos E 1 

The cartesian coordinates are then given by the relations 

(X P X \ 

= K z (n)K x (i)K z (u) \Y Py], (11) 
\0 0/ 

where Tv!^ and 7?. z describe rotations around the x- and z-axis, given by 

T^x(i) = | cosi —sin i | , Tl z {^) = | sinw cosw | . (12) 
cos ij 

It is then straightforward to show that the Hamiltonian (J3j) takes the form 

H = ~^ (13) 
and thus the equations of motion are given by 

A = M = 

A 3 

G = w = (14) 

K = tt = 0. 

Besides the actions A, G, K, the two angles u and f2 are also constants, which reflects the 
high degeneracy of the hydrogen atom. Equation (13) also gives a physical interpretation 
of A as a function of the energy. In quantum mechanics, A corresponds to the principal 
quantum number. 

Let us now return to the Zeeman effect. The Hamiltonian @ can be expressed in 
Delaunay variables as 

Hi = -^r 2 [l + cos 2 i + sin 2 i cos(2w + 2v)] , 

where r, sin 2 and v can be expressed in terms of Delaunay variables using (|6|), (||) and 
(|9|). The equations of motion take the form 

A = B 2 {A;H 1 } M = ^ + B 2 {M;H 1 } 

G = B 2 {G;H 1 } u= 5 2 {cj;#i} (16) 

K = ti = - + B 2 {Q;H 1 }, 
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where the Poisson bracket is defined by 

r, i _ df dg df dG df dg df dg df dg df dg 

U,9i dMdk dAdM dudG dGduj dVLdK dK dVt' 1 ' 

We discuss the computation of these Poisson brackets in Appendix [b| (see in particular 
Table i). 



To first order in B, (16) describes the Larmor precession of the ellipse. Since the 



Hamiltonian does not depend on f2, K = L z is a constant of the motion and (16) is in 
effect a two-degrees-of-freedom system. For small B, M is a fast variable, while A, G and 
uj are slow ones. The motion can thus be imagined as a fast motion of the electron along 
a slowly "breathing" and rotating ellipse. The dynamics can be visualized by a Poincare 
map, taking for instance a section at constant M, and plotting the value of G and oj at 



each intersection [DKN83]. 



2.2 Averaging 

To analyse the motion of ( |i~6| ) for small B, one can use the fact that M is the only fast 
variable, so that the dynamics of the slow variables will be essentially determined by the 
average effect of M during one period. 

With any given function / of the Delaunay variables, let us associate its average 

(f) A (G, K; u,n) = ±- ^ /(A, G, K; M, uj, Q) dM. (18) 
27T Jo 

The averaged Hamiltonian (H)\ generates the canonical equations 

A = M = ^3 + B 2 {M; (Hi) A } 

G = B 2 {G;(H 1 ) A } u= B 2 {u; <fTi) A } (19) 

k = o n = j + B 2 {n-,(H 1 ) A }. 

Since {Hi) a does not depend on M and £1, {H)\ is in effect a one-degree-of-freedom 
Hamiltonian, depending on K and A as on parameters. 

A standard result from averaging theory (see Appendix [A]) states that the equations 
( |j~9| ) are a good approximation of the equations (|i~6| ) , in the sense that 

• orbits of @ and (||) with the same initial condition differ by a term of order B 
during a time of order 1 / B 2 ; 

• to each nondegenerate equilibrium of (|l9|), there corresponds a periodic orbit of (li~6|), 
at a distance of order B 2 of the equilibrium, and which has the same stability if the 
equilibrium is hyperbolic. 

Averages over M can be computed quite easily by the formula 

1 f 2n dM 1 [ 2lx 

(/)a= 2^/ f{M) JE dE= 2^J f( M ( E ))(l-eco S E)dE. (20) 

Some useful averages are given in Table |l[ The averaged Hamiltonian can be written in 



the form |CDMW87| 



{Hi) A = ±A 4 (1 + cos 2 i)(l + fe 2 ) + |e 2 sin 2 icos2a; , (21) 
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</>A 


/ 


(/>A 


X 


-fae 


X 2 


a 2 (± + 2e 2 ) 


Y 





Y 2 


a 2 (i - ie 2 ) 


r 


a(l + ±e 2 ) 




a 2 (l + fe 2 ) 


z 


— fae sin i sin a; 


z 2 


a 2 sin 2 i [± + ±e 2 (3 - 5 cos 2a;)] 



Table 1. Some important quantities and their averages over the fast variable M. 

from which we deduce the relevant equations of motion 
G = 4£ 2 A 4 e 2 sinisin2w 

lb 



3(e 2 — 1) — 5 cos 2 % + 5(e 2 — sin 2 i) cos 2a; 



(22) 



The system can now be studied by analysing the orbits of ( p2[ ) or, equivalently, the level 
lines of the function ( |2T| ) in the (a;, G)-plane [DKN8S]. Doing this, one finds that the 
vector field has two elliptic equilibrium points located at 



2'T' 



G 2 = 



(23) 



which exist if \K\ < A/y/b. 

The vector field behaves in a singular way at the boundaries G = \K\ and G = A of 
the domain. These singularities have been explained by Coffey and coworkers [ CDMW87fl , 
who showed that the phase space has the topology of a sphere. The line G = A,ue[0, 2ir) 
has to be contracted into the north pole of the sphere, corresponding to e = and thus to 
circular motion; indeed, in this case the perihelion, and hence its argument ui, is undefined. 
The line G = \K\, uj £ [0, 2ir) has to be contracted into the south pole of the sphere, 
corresponding to i = and thus to equatorial motion; in that case, the sum + uj is 
sufficient to specify the position of the ellipse. 



To account for the spherical topology of phase space, [CDMW87] have introduced the 
variables 



£i = GAesini cos a; 
£2 = GAesini sin a; 
6 = G 2 -i(A 2 + ^ 2 ), 



(24) 



which belong to the sphere 



Using the fact that {H\)\ can be put into the form 

A 4 r 

{Hi) a = — 1 + cos z i + 3e z + 
lb L 



A 2 G 2 



(25) 



(26) 
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Figure 2 . Phase portraits of the averaged Hamiltonian ( po|) on the sphere £ 2 + £2 + £3 = 
(i(A 2 - A' 2 )) 2 , (a) for \K\ > A/VE and (b) for \K\ < A/y/E. The north pole corresponds 
to circular C-orbits, the south pole (not shown) to equatorial B-orbits, which are elliptic 
in both cases. In the second case, two additional periodic orbits (the Z-orbits) appear in 
the plane £1 = 0. 



and the Poisson brackets in Table || of Appendix [B|, they derive the equations of motion 



£1 = 


^6 


r Hi 

_A 2 G 2 


— 1 + e 2 + 5 cos 2 i 


6 = 




r Hi 
-A 2 G 2 


-4(l-e 2 ) 




£3 = 


5A 2 








AG 







(27) 



In these variables, the poles have become equilibrium points around which the flow is 
nonsingular, so that their stability can be easily determined. Depending on the relative 
value of the constants of motion A and if, there are two qualitatively different phase 
portraits (Fig. ||): 

1. If I If I > A/y/b, both poles are elliptic, and there are no other equilibrium points. All 
other orbits rotate around the sphere with Co < (Fig. ||a). 

2. If I if I < A/v5, the south pole is still elliptic, but the north pole has become hyperbolic, 
and the new elliptic equilibria (|23|) have appeared in a pitchfork bifurcation. Two 
homoclinic orbits of the north pole separate orbits rotating around each of the three 
elliptic equilibria (Fig. |2|t 



The averaging theorem shows that to each of the four possible equilibrium points 
of (|27D, there corresponds a periodic orbit of the exact system ( |l~6| ) (see also Appendix 
|B|). For further reference, let us call B-orbits the equatorial orbits (which lie in the 
plane perpendicular to B), C-orbits the circular ones, and Z-orbits those corresponding 
to the nontrivial equilibrium (p3|). One further expects that the periodic orbits of ( |27| ) 
approximate either quasiperiodic KAM-type orbits of (|l~6|), or "soft" chaotic components 
associated with resonances. More prominent chaotic components are expected near the 
homoclinic orbits of the north pole. 

In other works [HRW89], orbits are sometimes represented in cylindric coordinates 
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(p, (f), z). They can be deduced from Delaunay variables by the relations 

z = rsini sin((jj + v) 



p = ry 1 — sin i sin (u + v). 



(28) 



The different periodic orbits can thus be parametrized either by the true anomaly v or by 
the eccentric anomaly E as 



C-orbits z = A\/ A 2 — K 2 sin v p = A\/ A 2 cos 2 v + K 2 sin 2 v 

B-orbits z = p = A(A - \/A 2 - K 2 cos £) (29) 



Z-orbits z = A 2 sin i (cos E — e) p = A 2 y (1 — e cos E) 2 — sin 2 « (cos E — e) 2 , 

where sin 2 i = 1 — |iT|/\/5A and e 2 = 1 — \/5|if|/A in the last case. The C-, B- and 
Z-orbits are labelled, respectively, C, I\ and Iqo in 



3 The Stark effect 

We consider now the Hamiltonian (|T|) in the case B = 0. If we choose the z-axis along the 
electric field F, it can be written as 

H = -p 2 - - + Fz. (30) 
2 r 

Besides the energy and the z-component of the angular momentum, this system has a 
third constant of the motion and is thus integrable. Indeed, as shown by |R63 1 , the 
generalization of the Runge-Lenz vector (Q) , 



satisfies the equation of motion 



and thus C 7 = C • F is constant. 



A — ^ (r A F) A r (31) 



C = ^LAF, (32) 



We will start, in Section 3.1, by analysing the system in Delaunay variables, in particu- 



lar in its averaged form, in order to get a feeling for the geometry of the orbits. In Section 



3.2| , we present another description of the system, based on parabolic variables, which 
allows to construct action-angle variables taking the constant C z into account. Though 
they are better suited for perturbation theory, these action-angle variables have a less 
obvious geometric interpretation than Delaunay variables. This is why we establish the 
transformation formulas between both sets of variables in Section |3.3| , in the limit F — > 0. 

3.1 Delaunay variables 

In Delaunay variables, the Hamiltonian (30) takes the form 



H = ~ 2X 2 " + ™ (A ' G ' K] M > W) ' (33) 
H2 = r sini sin(a; + v), 
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7T 2tt 

Figure 3. Phase portraits of the averaged Hamiltonian ( |36| ) (a) in the (w,G)-plane and 
(b) on the sphere £ 2 + £f + £f = (^(A 2 — AT 2 )) 2 , seen from below the south pole. The 
elliptic points at u> = ? , ^ correspond to periodic S-orbits of the Stark Hamiltonian. 



and the equations of motion have the structure 
A = F{A;H 2 } 
G = F{G;H 2 } 

k = o 

The constants of motion are H, K and 
C 2 = — esini sinc<j 



M = ^ + F{M-H 2 } 
uj= F{u;H 2 } 

n= F{n-,H 2 }. 



(34) 



—Ft 1 [l — sin 2 i sin 2 (w + u )] . 



(35) 



In order to understand the geometry of the orbits for small F, we may analyse the averaged 
Hamiltonian 



(H)a 



1 



2A 2 



— Fae sin i sinw. 
2 



(36) 



The relevant equations of motion of this one-degree-of-freedom system are 

G = |i ? A 2 esinz cosw, 



, G A -K 2 K 2 

w = i F ^3 ■ ■ 

G^esi^ 



(37) 



sinw. 



We observe the existence of a pair of elliptic stationary points at u> = |, ^ and G 2 = |if |A, 
which implies e = sini = ^1 — |i^|/A (Fig. |3]a). When F = 0, the constant of motion 
C 2 reaches its extremal values ±(1 — K/|A|) on these points. We will call S-orbits the 
associated periodic orbits of the Stark Hamiltonian. In order to analyse the motion at the 
boundaries of phase space, we use again the variables (|24|). Since the averaged Hamiltonian 
can be written as 



3 A 
2G 



6, 



(38) 
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Figure 4. Same as Fig. [| but for K = 0. 



they evolve according to 



£1 

e' 3 



l^[-^-(A 2 + ^ 2 + 2e 3 )e3] 



3FA&. 



1^ 7^66 



(39) 



This shows that there are no other equilibrium points than the two elliptic ones, since the 
flow is nonsingular at the poles (Fig. ||b). 

Besides the periodic orbits corresponding to the elliptic equilibrium points of (|37|), the 
system displays quasiperiodic orbits for which both uj and G oscillate. There is a particular 
orbit following the meridian £2 = of the sphere, for which the Kepler ellipse oscillates 
between a circular and an equatorial one, while its major axis is always perpendicular to 
F. Note that the case K = is special, since the elliptic points merge at the south pole. 
All orbits then go through the south pole, where the eccentricity vanishes, which means 
that the electron approaches the nucleus arbitrarily closely (Fig. 



3.2 Parabolic and electric action— angle variables 

The separability of the Stark Hamiltonian in parabolic variables was already known to Max 



Born [B27] from the earlier works of P.S. Epstein [E16] and K. Schwarzschild [SchwlC]. 
Parabolic variables (P^, P^, P v ; £, rj, ip) are defined by 



x = £77 cos ip 
y = £7/ sin ip 



Pa 



Py 



Pz 



— f- — — cos ip - — sm ip 



£2 + ^2 

& - vP, 

£2 + ^2 



sm p> + — — cos ip 



(40) 



where P v = L z = K is the z-component of the angular momentum. 



In these variables, the Hamiltonian takes the form 



H 



1 



rl 



C + rj 2 



K 2 



-Pt + — + Ft 4 + -P z + ^~ 



K 2 



Fif 



(41) 



11 




FIGURE 5. Level lines of the constants of motion a± and ol-i for H = 1, K = | and 
F = 0.05. In contrast to £, r\ may have an unbounded motion. However, the saddle point 
is located quite far away from the origin, at rj 2 = 0(F^ 1 ). 



There are four constants of motion H, cti, ot^, K, related by 

1 rj K 2 



Oil 



a.2 



2 p i + ^-^ H+ F ? 



1 



2' r > 
Oil + «2 = 2. 



K 2 
2r/ 2 



rj 2 H- Frf 



(42) 



If one scales time by a factor £ 2 + rf , ( |4l| ) is seen to describe the motion of two decoupled 
"oscillators". In fact, the constant ai is of the form \P 2 + Vi(£), where Vi(£) — > oo m 
both limits £ —* and £ — » oo when F > 0. Thus the motion of £ is always bounded. By 
contrast, c*2 = + ^(ry), where 1^(77) — ► —00 for T] —> co. Thus the level sets of 02 
may be unbounded for large values of r? 2 or «2 (of order i 7-1 ), corresponding to ionization 
(Fig. H). The following discussion is limited to the bounded motion of 77, which exists for 
small F. In that case, action-angle variables can be constructed. 

Action variables related to the constants «i and ai are defined by 



■h 



2vr 



'2H + 



J,? " 27T 



2a 2 K 2 



(43) 



where the integrals are over bounded level sets of a\ and ol^ respectively. In the case 
F = 0, they can be computed explicitly [B27]: 



J, 



1 



-K + 



Oil 



^2H 



■J f! 



-K + 



a- 2 



From this we deduce that 
H = - 



«1,2 



^2H 

2J(, V + K 



(44) 



(45) 



2(Jf + J„ + i^) 2 ' Jt + J v +K 7 

which shows, by comparison with (|T3|), that Jg + J v + ivT = A is nothing but the first 



(46) 



Delaunay action. Moreover, one finds that the constant of motion C z reduces to 



1 — OL2 = Oil — 1 



A 
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which suggests to introduce the action J e = J v — J^. In quantum mechanics, J e corresponds 
to the electric quantum number. For further reference, let us call (A,J e ,K) the electric 
action variables. As we have seen that \A Z \ ^ 1 — |-fC|/A for F = 0, they vary on the 
square domain 

-(A- \K\) ^ J e ^ A- \K\. (47) 

For F > 0, we need to know the expression of the Hamiltonian in action variables. 
This can be done perturbatively [DK83] with the result 

H(A, J e ,K) = _ 3FAJ e - ^F 2 A 4 (17A 2 - 3J e 2 - 9K 2 ) + 0(F 3 ). (48) 

2A Z 4 

The associated canonical equations are 

A = w A = ^ - 3FJ e - §F 2 A 3 (17A 2 - 2J 2 - 6K 2 ) + 0(F 3 ) 



j c = We = -3FA+ §F 2 A 4 J e + 0{F A ) (49) 

K = w K = %F 2 A 4 K + 0(F 3 ), 

where w\, w e and wk (the electric angle variables) are conjugated to A, J c and K re- 
spectively. In the case F = 0, these equations are equivalent to the equations (14) in 
Delaunay variables. The electric field suppresses the degeneracies of the Kepler problem, 
and introduces different time scales for the various angles. 

In order to compute the electric angle variables, we need to parametrize the level curves 
of Oi\ and 02- For F = 0, using (|45|) the first equation of (H3) can be written as 



, 2 , (e 2 -exf -a) 

A 2 £ 

where ^ = a\ ± 61 are extremal values of £ 2 with 



+ - — — — = °> ( 5 °) 



ai = A(2J ? + K), h = 2A^J^ + K). (51) 

Similarly, ?? 2 can vary between 02 — ^2 and 02+^2- In terms of the electric actions (A, J e , if), 
these limits are given by 



ai = A(A - J e ) 61 = A V / (A - J e ) 2 - if 2 

a 2 = A(A + J e ) 62 = A^(A + J c ) 2 -K 2 . (52) 

Relation (|50|) is the equation of an ellipse in the (£ 2 , A£i-g)-plane, which suggests to 
parametrize the level sets by 



£ = -^/ai - 61 cos V> f? = V a 2-b2 cos x 

D Mini/* □ 62 sin x , K „, 

P « = ~W~ ^ = ~An~' (53) 

where ip and x are auxiliary angles, playing a similar role as the eccentric anomaly E in 
the case of Delaunay variables. If we introduce the action 

S = y ? P 5 de' + P v dr]' + Kip, (54) 
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the angles conjugated to (Jg, J V ,K) are given by the formulas 



dS 



w„ 



dS 



dS 
dK' 



(55) 



and the angles conjugated to (A, J C ,K) are then obtained by the linear transformation 

»A = l(ui£ + w v ), w e = ^(w v -w^), WK = W ip -W\. (56) 



With the parametrization (|53[), the derivatives (p5j) take a simple form [B27|, and the final 
result is 



V> + X 
2 

x- i> 



1 



2A 2 



(6i sin ip + 62 sin x) 



(57) 



W K = ¥ ~ Pl(V0 - P2(X), 



where 



Pji e ) = — 



69' 



a j 



bj cos 9' ' 



1,2. 



(58) 



Using the fact that a\ — b\ = K 2 A 2 , this expression can be written as 
Pj(9) = sign(ET) atan 4 / - ' / tan • 



aj - bj 



for — 7T ^ 6 ^ 7T, 



(59) 



and can be continued to arbitrary 6 by the rule + k2ir) = Pj(9) + kir for A; G Z . This 
implies that for all 9, we have 



cos Pj(9) 



cos - . 

bi cos 9 2 



J aj + bj 9 
sin pj (9) = sign if— ^= „ sm o- 
aj — bj cos 9 £ 



(60) 



Note in particular that the transformation w e 1— * w e + ir, w\ 1— * + ir corresponds to 
keeping ip fixed and increasing \ by 2-7T. Hence it leaves all parabolic variables fixed, 
except cp which is increased by it. In other words, this transformation describes a rotation 
of angle ir around the z-axis. 

The relations ( |57D are valid for F = 0. Higher order expressions in F of the angle 
variables can be computed perturbatively, see Appendix |C|. 

3.3 Correspondence between electric action— angle and 
Delaunay variables 

We now establish transformation formulas between electric action-angle variables and 
Delaunay variables in the case F = 0. In the two-dimensional case, corresponding to 
K = 0, this issue has been addressed in [DW91], in connection with Lissajous variables. 

The most important relations can be obtained by averaging over the fast variable u>a- 
In view of d57j), the averaging operation can be written as 



</(V>,X)>A 



1 

2^ 



2k 



1 



61 cos ip + 62 cos(?/> + 2w e 



2A 2 



dip. (61) 
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/ 


</>A 


f 


(/)a 


cos ip 
simp 
cos 2 ip 
cos ip cos x 
sin ip sin x 


-^-^cos2w e 

jfc sin 2w e 
1 

2 

i cos 2ui e 
i cos 2w e 


cosx 
sinx 
cos 2 x 
cos ip sin x 
sin ^ cos x 


-^cos2 Wc -^ 

sin2w e 
l 

2 

— ^ sin27i; e 
i sin 2w e 



Table 2. The averages over the fast variable wa of some important quantities. Using 
( |53| ) they can be used to compute averages of some polynomials in £ 2 and rj 2 . 



A few useful averages are given in Table We can thus easily compute the averages 
1 



fAJe 



1 



<r>A= §<£ 2 + »7 2 >a 



1 + 



1 



8A 4 



(b\ + b 2 2 + 2hb 2 cos 2w ( , 



(62) 



Comparison with the corresponding averages over M (Table |) gives us the relations 

J e = Ae sin i sina;, 



1 



4A 4 



(6 2 + b\ + 26i6 2 cos2w e ). 



(63) 
(64) 



This last relation allows to compute G and i. It has a nice geometric interpretation: the 
maximal value of e 2 (hence the minimal value of G) is attained for w e = and 7r, while 
the minimal value of e 2 is reached for w e = ^ and According to Fig. ||a, these values 
of w e correspond to w = | sign J e . From the fact that w e < for F > 0, we infer that 



sign(cosa;) = — sign(sin 2w 



(65) 

The relations (|6"3|), (|64|) and ( |65| ) determine the transformation (G,u) ^ (J e ,u> e ) up to a 
phase 7r of u! e , which depends on S7. 

It remains to establish relations between the angles (wa, wk) and (M, fi). Since wa = 
M for F = 0, the difference u>a — M does not depend on M. From ( |53D we deduce that 
in electric action-angle variables, 



A 2 



1 



1 



2A 2 



?i cos ip + 62 cos x) 



A 2 



<9u>A cfotv 



<9^ 5x 



Comparison with (||) yields the equalities 

1 



and, using wa = 



_d_ 

dE 



ecosE 



dM d 



2A 2 



(61 cos ip + b 2 cosx), 



dE dM \ dip 



rdwA 9wa\ d d d 
dx ' dwA dip dx 



(66) 



(67) 



(68) 
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Applied to (|6^) , this also gives 

esvaE 



2A 2 



?i sin ip + b 2 sin x) • 



(69) 



In particular, when M = we have E = and thus by ( |57| ) j«a = Inserting the 



relations V = — w e and x = w\ + w e into (^7|) and (^) , we can solve for cos w\ and 
sinwA with the result, for general values of M, 



cos(wa — M) 



2A 2 e 



cos w c , 



sin(?z;A — M) 



2A 2 t 



h + b 2 

Relation (|68|) also implies that 

ip — E = u>a — w e — M, x — E = w\ + w e — M. 



sin w e . 



(70) 



(71) 



Determining wk is a bit more delicate. We will use the fact that in Delaunay variables, 
the x-component of the angular momentum is 



L x = yjG 2 - if 2 sin ft, 
while in parabolic variables, we have from ( |40| ) 



1 IK 

L x = gfa-Ffe - C^)sin^ - ~ *7 2 )a»¥>- 



(72) 



(73) 



Being independent of M, L x also has to be independent of wa. We may thus evaluate 
(|73|) in the case tf) = 0, x = 2w e . Then (|5"3|) reduces to 



£ = yoi — 61, ?? = y / a 2 -b 2 cos2w e , ^ = 0, 
and ([[?]) implies = wk + P2, where 



Pn 



b 2 sin 2w e 
Ai] 



cos p 2 



~ b 2 , . , „s VQ2 + 62 . 

cos iUe , sm /32 = sign( K ) sin w e . 



>1 



Inserting (J7J) and (|75|) into (73), we obtain 

L x = L\ cos wk + £2 sin u;^, 

where 

Li 1 



(74) 



(75) 



(76) 



cos w e 



2A^ 2 



2sign(i^)^V«2 + b 2 b 2 sin 2 w e + KA{£ 2 - r) 2 )^ a 2 - b 2 



\ K 02 - b 2 

2 V «i - h 2Ari 2 ^a 2 - b 2 



1 



KA(a 2 -b 2 ) + 2 sign(K)b 2 J a\ - b 2 , sin 2 w, 



-K 



a 2 - b 2 
a\ — b\ 



n 



a 2 - b 2 



(77) 



where we have used the relation \/ a 2 , — b\ = |-K"|A. The term L 2 can be evaluated in 
a similar way. We obtain the following, relatively compact expression of L x in electric 
action-angle variables: 



L x = Yi(A, J e , K) cos w e cos wk + Y 2 (A, J e , K) sin w e sin wk, 



(78) 
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where we have introduced the notations 



y 2 (A, j c ,k) 









V 


2 




1, , 


- 2 \K\ 





ia 2 



a\ — b\ 



CL2 



a 2 + b 2 



a i 



(79) 



ai-bi V a 2 + &2 
Comparison with ( |72|) gives the desired relation between and O: 

VC^p sm (^ _ = Yl (A, J e , if) COS W e 

Y 2 (A, J e , K) sinu; e . 



y/G 2 -K 2 cos (SI 



(80) 



We point out that despite the absolute value, the expressions ([79]) are regular at K 
and admit the Taylor series 



Y 2 (A,J C ,K) 



-K- 



Je, 



+ K- : 



A 2 X 



2(A 2 - J|) 5 /2 

A 2 



+ 0(K 5 



(81) 



2 (A 2 - J e 2 ) 3 / 2 



+ 0{K 4 



4 The crossed-fields problem 



We now consider the full crossed-fields Hamiltonian ([l]), in the case < B, F <^ 1. First, 
we have to choose a system of coordinates. Both sets of action-angle variables that we 
have used so far are defined with respect to a privileged direction (the z-axis). Depending 
on the regime we consider, it will be most convenient to choose this direction along the 
electric or along the magnetic field. In the first case, the Hamiltonian takes the form 



H 



-p 2 -- + Fz + -BL X + -B 2 (y 2 + z 2 ). 
2 F r 2 8 v ' 



(82) 



To account for the second case, we also introduce coordinates (x',y',z') = (z,y,x). In 



Section 4.1, we consider the case B <C F, which is a small perturbation of the integrable 
Stark effect, and thus particularly well suited to perturbation theory. The case F <C B 



is considered in Section 4.2. The orbits contained in the plane perpendicular to B exist 
for all values of the fields. We analyse them in Section |4.3| . Other periodic orbits and the 



general structure of phase space are discussed in Section 4.4 



4.1 The case B < F 

When the magnetic field acts as a small perturbation of the Stark Hamiltonian, it is best 



to use the electric action-angle variables introduced in Section ^2. The Hamiltonian can 
be written as 



H = H (A, J e , K; F) + BH^A, J c , K; w e , w K ; F) 

+ B 2 H 2 (A, J e , K; w\, w e ,w K ; F), 



(83) 



where Hn 



l 



3FAJ e + F 2 h 2 (A, J e , K ; F) is the Stark Hamiltonian (^8|), Hi = ±L X 



has been computed in d78|), and H 2 = |(y 2 + z 2 ). The equations of motion thus have the 
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structure 

A = 0(B 2 ) w A = -^-3FJ e + 0(F 2 ) + 0(B) 

A J 

j e = 0(B) w e = — 3FA + 0(F 2 ) + O(B) (84) 

K = 0(B) w K = 0(F 2 ) + 0(B) 

We can again average over the fast variable w A , using the rule @). Since L x does not 
depend on w A , Hi is already in averaged form. The equations of the averaged system are 
thus given by 

X= B{J C ;H 1 } + B 2 {J C ;(H 2 ) A } 

K= B{K;H 1 } + B 2 {K;(H 2 ) A } 

(85) 

w e = -3FA + F 2 {w e ; h 2 } + B{w c ; H x ) + B 2 {w e ; (H 2 ) A } 
w K = F 2 {w K ; h 2 } + B{w K ; Hx} + B 2 {w K ; (H 2 ) A } 

Since we assume that B <C F, w e evolves on a faster time scale than wk- We may thus 
further approximate the dynamics by averaging over w e , that is, we define the double 
average 

((f))A,j e (K;w K ) = J J f(A,J e ,K;w A ,w e ,w K )dw A dw e . (86) 

The doubly-averaged Hamiltonian {{H)) At j c has one degree of freedom and is thus inte- 
grable. Let us now compute various averages, at lowest order in F. In order to compute 
(H 2 ) A , we need to evaluate the average of 

2,2 ,-22-2 , 1//-4 , 4 r,<^2 2\ / ^\ 

y +z = C 7? sin (p + 3 (£ + rj - 2f 77 ). (87) 

Using the expressions in Table [2|, the second term is easily averaged. To average the first 
term, we use the fact that with (|57|) and (|60[) , ^Tjsin^ can be written as a polynomial 



in vVl i ^l) V fl 2 i ^2) and sines and cosines of iwr-, $ and ^. The final result after 
simplification is 

(y 2 + z 2 ) A = ±A 2 (2A 2 + 3J 2 - K 2 ) - |b x 6 2 cos 2w x 

+ 3 6i& 2 cos2w c - ±A 2 (A 2 - J 2 - ET 2 )cos2u; c cos2t<;K (88) 
- A 2 J e Ksin2u; e sm2w K + 0(F). 

Computation of the doubly-averaged Hamiltonian is now easy. From (|T^) and (|88|) we 
obtain 



((L x ))a,j c = O(F) 

^A (2^ i . , , j , 



((y 2 + z 2 )) AJc = ±A 2 (2A 2 + 3J 2 - K 2 ) - |6 1 6 2 cos2^ + 0(F). 



In fact, equation (122) in Appendix |C| shows that also at order F (and probably at all 
higher orders), the transformation w e i— > w e + 7r, wa ► + ^ describes a rotation of 
angle tt around the z-axis. Since this rotation changes L x into — L x , we conclude that 
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( a ) k 3 = K ( b ) k 3 ( c ) 




Figure 6. Phase portraits of the doubly- averaged syste m (J9C| ) on the sphere 
(A - | J 1) 2 , (a) for < B < B 2 , where B 2 is given by (J95|) , (b) for B 2 < B < B u given 
by (|93|), and (c) for B > B\. The poles (the intersection of the sphere with the K3-axis) 
correspond to the S-orbits present in the Stark effect. The intersections of the sphere with 
the K\- and K2-axis correspond respectively to the BF-orbits, lying in the plane defined by 
B and F, and the B-orbits, lying in the plane perpendicular to B. 



((L x ))\j a = 0(F 2 ). Discarding irrelevant constant terms, the doubly-averaged Hamilto- 
nian can be written as 

((H))a,j c = |F 2 A 4 K 2 - ±B 2 [A 2 K 2 + |6 1 6 2 cos2^] + 0(F 3 , BF 2 , B 2 F), (90) 

where b\ and 62 are given in ( |52[ ) and the remainder denotes a sum of terms of order F 3 , 
BF 2 and B 2 F. Besides A and J e , © is a third adiabatic invariant of the crossed-fields 
Hamiltonian in the case B <C F <C 1. Up to the remainders, the equations of motion have 
the form 



3 A 2 1 (91) 

l-^(A 2 + J 2 -^ 2 )cos2^ . 



3 2 

K = -—B~bib 2 s'm2wK 
lb 

w K = -F 2 K i K - -B 2 A 2 K 
2 8 

We should also examine the topology of phase space. From (|47|), we deduce that in the 
limit F — > 0, K varies between —(A — | J e \) and A — | J e \. If, say, J c ^ and K = A — J e , 
we have b\ = and the angle tp is undefined. Since pi(tp) = 5 in this case, we obtain 
from (|5^) and (p7|) that the quantities w\ + it?^ and w\ + w e are sufficient to determine 
the state of the system completely. We conclude that in the averaged phase space, the 
variables wk and w e are irrelevant when \K\ = A — | J e |, and thus we have again a spherical 
topology. The sphere can be parametrized by 



(«i,«2,« 3 ) = (^(A - I J e |) 2 - K 2 coswk, V(A - I Jc|) 2 - K 2 smwK, K). (92) 

We now analyse the structure of phase space for increasing B. When B = 0, the orbits 
of ( |9l|) follow the parallels of the sphere. The poles and all points of the equator K = 
are fixed points. For slightly positive B, a resonance of order 2 is created: only the points 
wk = 0, f ,7r, of the equator remain fixed. A straightforward stability analysis shows 
that the points (wk,K) = (0,0) and (tt, 0) are always elliptic, while the points (§,0) and 
(^f ,0) are hyperbolic for 



/A 2 — T 2 

B<Bl = 6V2M.L. (93) 
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and elliptic for B ^ B\. If J e 7^ 0, there is another pair of equilibrium points, located at 

7T 37T 
2 ' 2 



wr = § j tt and = , given by the condition 



B - 1 2v/[(A- J.) 2 -if, 2 ][(A + J c )2-ir2] 

These orbits are created in a pitchfork bifurcation at the poles at B = (Fig. ^a), move 
to the equator as B increases (Fig. ^b), and disappear in another pitchfork bifurcation at 
B = B\ (Fig. ||c). It is easy to see that there must be a global bifurcation involving a 
saddle connection between these values. It is given by the condition 



((H)) A;Je (Z,0) = ((H)) AiJe (w K ,A-\J e \) B = B 2 = 6V2FAJ ^-^ . (95) 

At B = B2, the poles of the sphere are connected with the points (§, 0) and (^f-, 0) on the 
equator by heteroclinic orbits. In the case J e = 0, all these bifurcations collapse. 

Our analysis of the doubly-averaged Hamiltonian (|9^) has thus revealed a rather rich 
structure of phase space. We point out that only the first case, depicted in Fig. |6]a, 
is compatible with the hypothesis B <C F, which is necessary for the doubly-averaged 
system to be a reliable approximation. We will see in the next sections, however, that the 
picture given in Fig. ^ also contains some truth in the other parameter ranges (see Fig. 

n. 

We conclude that for B <C F, the structure of phase space is determined by four types 
of orbits: 

1. The poles of the sphere ( p2| ) correspond to the fixed points of the averaged Stark 
Hamiltonian, and thus to the periodic S-orbits of the original (unaveraged) Hamilto- 
nian (|83|). They are unstable unless J e = 0. 



2. The points K = 0, wk = § or t?, which are hyperbolic for small B, correspond to 

= \ or ^ in Delaunay variables, and hence describe the B-orbits, which lie in 
the plane perpendicular to B. In the singly-averaged system, they appear as periodic 
orbits, following a curve of constant J e . Since K = 0, all these curves agglomerate 
at the points G = and uj = ^ or ^ (see Fig. ||). In the original system, the 
orbits can be interpreted as a fast rotation along a slowly "breathing" Kepler ellipse, 
reaching periodically the eccentricity e = 1, where the electron approaches the nucleus 
indefinitely closely. 

3. The points K = 0, wk = or tt, which are elliptic, correspond to O, = or it, and 
thus describe orbits in the plane of F and B. As in the previous case, they evolve on 
a level curve of J e containing a point with zero angular momentum. We will call them 
the BF-orbits. 

4. Finally, the four points (if, i-K*) describe more complicated orbits, which are stable 
and provide a connection between S-orbits and B-orbits. For small B, they are close 
to the S-orbits. Let us thus call them SB-orbits. 

4.2 The case F < B 

When perturbing the Zeeman effect, it seems more appropriate to use coordinates (x' , y' , z') 
in which the magnetic field is vertical. The Hamiltonian takes the form 

H = \p 2 -- + \bL z , + l -B\x 12 + y' 2 ) + Fx'. (96) 
2 r 2 8 
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We will denote by (A', G', K'; M', u>', ft') the associated Delaunay variables. The trans- 
formation between these Delaunay variables and those defined with respect to (x, y, z) = 
(z', y', x') can be derived by expressing L and z in both sets of variables. The result is 

A = A' G = G' 

K = a/G" 2 - K' 2 sin ft' sin 2 i = cos 2 ft' + cos 2 % sin 2 ft' 

\J G' 2 — K' 2 cos ft' sin uj' cos ft' + cos i' cos u>' sin ft' 
cos ft = — ^^^^^^=^^^^^^= cosu; = (97) 

VG' 2 cos 2 ft' + K' 2 sin 2 ft' sm * 

K' cos oj' cos ft' — cos i' sin <J sin ft' 
sin ft = — ^^^^^^^=^^^^^^^ sin u = . 

VG" 2 cos 2 ft' + K' 2 sin 2 ft' sm * 
The averaged Hamiltonian over M' takes the form 

{H)A = ~ 2X 2 + I K + ^ 2 ^i)a( G '' + F(H2)a(G', K>; J, ft') 

(H^ A = ± [(1 + cos 2 0(1 + |e' 2 ) + |e' 2 sin 2 i' cos 2w'] ( 98 ) 
(-^2) A = — § A 2 e' (cos u/ cos ft' — sin uJ sin ft' cos i') . 
It generates the equations of motion 

K'= F{K';(H 2 ) A } n' = ^ + B 2 {n , ;(H 1 ) A } + F{n';(H 2 ) A } 

G' = B 2 {G';(H 1 ) A } + F{G';(H 2 ) A } u' = B 2 {uj'; (Hi) A } + F{uj'; (H 2 ) A }. (99) 

We observe that for F <C B, ft' evolves on a faster time scale that a/, and the dynamics 
can be further approximated by averaging the Hamiltonian over ft'. Note, however, that 
the average of {H 2 ) A over ft' vanishes, while {H\) A does not depend on ft'. The twice 
averaged Hamiltonian is thus strictly equivalent to the averaged Zeeman Hamiltonian ( |2"T| ) . 
This means that to lowest order in perturbation theory, the electric field does not influence 
the phase portrait of the Zeeman effect, which has one of the two behaviours indicated in 
Fig. [2]. The phase space is thus organized around three types of periodic orbits: 

1. a stable B-orbit at K' = G", located in the plane perpendicular to B; 

2. a circular C-orbit, stable for large K' and unstable for small K'; 

3. and a pair of stable Z-orbits, existing only for small K' , with a major axis perpendicular 
to the line of nodes. 

These orbits will experience deformations of mag nitude 0(F 2 ) when F > 0. 

4.3 The B-orbits 

Let us now examine the most important periodic orbits for general (small) values of F and 
B. A special role is played by the orbits in the plane perpendicular to B, which we called 
B-orbits. Their existence can already be seen on the original Hamiltonian (p^), which 
leaves the plane x = invariant. The Hamiltonian restricted to this plane has two degrees 
of freedom, its averaged version will thus have one degree of freedom. Hence all orbits 
starting in that plane will look periodic in the averaged system, but may correspond to 
quasiperiodic or soft chaotic components of the original Hamiltonian. 

In the two previous sections, we have found that the B-orbits are hyperbolic in the limit 
B — > and elliptic in the limit F — ► 0. Hence there must be at least one bifurcation value 
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(a) 
A 







(b) 
A 



2?r 







2?r 



Figure 7. Orbits of the averaged Hamiltonian (IOC) restricted to the plane perpendicular 
to B. (a) shows a case with B < F, and (b) shows a case with F < B. There are two 
types of orbits, those which cross the lines w c = and tt, and those which oscillate around 
w c = 5 and t?. The second type will not appear on a Poincare section at w c = 0. 



between these limits. The twice averaged Hamiltonian in electric action-angle variables 
fl90|) suggested that this transition should be given by the condition (93), which is, however, 



not in the range where (|90|) can be expected to be a good approximation. We will now 
examine this question in more detail with the once averaged Hamiltonian. 



Using the expansion (81) of L x for small K, it is easy to compute the equations of 
motion (|8^) of the singly-averaged Hamiltonian in electric action-angle variables. One 
can then check that both wk and K vanish for K = and cos Wif = 0, which confirms 
the invariance of the subspace of B-orbits. The motion in this subspace is determined by 
the one-degree-of-freedom Hamiltonian 

(H) A (J e , 0; w e , f ) = - 3FAJ C + f F 2 A 4 J c 2 + 0(F 3 ) 

- \B^A 2 - J 2 sin u; e + ^B 2 A 2 [A 2 + J 2 + (A 2 - J 2 ) cos 2w e ] 
+ 0(BF). (100) 

Up to the remainders, its equations of motion are given by 

j e = -B^A 2 - J| cos u; e + — B 2 A 2 (A 2 - J 2 )sm2w e 

3 1 J 3 ( 101 ) 

w e = -3FA + -F 2 A 4 J e + -B , ° = sin w e + -B 2 A 2 J e sin 2 w e . 
2 2 ^A 2 - J e 2 8 

This system admits two elliptic equilibria located at w e = \ and 4f and 



±A 



/ B V 
1+ (6?a) 



-1/2 



(102) 



These points move from the boundaries | J e \ = A to J c = as % increases from to oo. 
The orbits of (|101| ) are shown in Fig. ^. They are of two types: 

1. the orbits which cross the lines w e = and n, and for which w e is monotonous; 

2. the orbits which oscillate around w e = \ and 4^ without reaching w e = or tt. 

There is no separatrix between both types of orbits, the flow being regular when expressed 
in the right variables (the boundaries \J e \ = A should again be contracted into the poles 
of a sphere, since they correspond to the periodic S-orbits of the Stark effect). Still, there 
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Figure 8. The two types of equatorial orbits: (a) singular orbits, which periodically reach 
an eccentricity e = 1, and (b), regular orbits which are bounded away from e = 1. The 
electric field is vertical, while the magnetic field points out of the plane. One can identify 
the fast motion along a Kepler ellipse, the parameters of which evolve approximately on 
a level curve of Fig. || in the regular case, and of Fig. [| in the singular case. These orbits 
are not periodic in general, as can be seen on the pictures. 



is an important qualitative difference between both types of orbits. Indeed, by (p4[), the 
values w e = 0, 7r are the only ones that lead to an eccentricity e = 1 when K = 0, and thus 
imply a close encounter with the nucleus. Hence, the first type of orbits will contain cusps 
(singular orbits, Fig. |8|a), and the second will not (regular orbits, Fig. |8]b). The boundary 
between both types of orbits, when starting on the line w e = ^, is given by the condition 



<fr>A(A,0;f,f) = (fr>A(Je,0;f,f) 



B ~ 6FA 



A- J e 

A + J c 



(103) 



We will now examine the stability of the 5-orbits in the 4-dimensional phase space of 
the averaged Hamiltonian. To do this, we first need to compute the derivatives 



du)K 



dwK 
dwK 



5.0 2 



dK 



dK 



B — ; ^ C cos w c H — B 2 A 2 J e sin 2w c 
VA 2 - Jf 4 

9 rfA4 B A 2 

=-F A +— — -r—prsmw e 

o 2 2 (A 2 - J 2 ) 3 / 2 



dK 
~dK 



5,o 



A 2 



1r 

16 A 2 -J 2 



(104) 



[5 A 2 + J 2 + (3A 2 - J 2 ) cos 2w c 



8wk 



- 

2 ' u 



-B^A 2 - Jf sin u; e + -5 2 A 2 (A 2 - J 2 ) [3 + 2 cos 2w c 



In order to determine the stability of a B-orbit, we have to find the multipliers of the 
variational equation 



z = A(w e , J e )z, 



(105) 



integrated over a solution of the system ( |101 ), where z T = (wk — §>0)> an d A is the 
matrix with entries given by (104). 

Two limiting cases can be studied analytically. When B — > 0, A can be replaced by 
its average over the fast variable w e , and we obtain again the condition (|93|), which tells 
us that only orbits with A 2 — J 2 = 0(B 2 / F 2 ) are elliptic. The other case is that of the 
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1 2 0.0 0.5 1.0 

FIGURE 9. Stability diagrams of the equatorial orbit, obtained by computing numerically 
the multipliers of equation ( |105| ), starting with w e = § and different values of J e , F and B. 
(a) shows the case J c — (b) the case F = 0.1. Shaded regions indicate unstable orbits. 
The curves are: I the location of the stable orbit (102), II the boundary (|103[) 



between 

regular and singular orbits (orbits with large B, large J c or small F are regular), and III 
the stability boundary (B3|) of the doubly-averaged system. 



orbits (|102|) , for which both w e and J e are constant. In this case, we find det^4 > 0, which 
means that the orbit is elliptic. 

Fig. |9] shows numerically computed stability diagrams for general parameter values (one 
should note that equation ( |105| ) is much easier to treat numerically than the full equations 
of motion). As expected, orbits are unstable in a region compatible with condition (|93|) 
and stable for B 3> F. There appears to be a second zone of instability for intermediate 
values of -p. However, this region corresponds to rather large values of B, for which the 
averaged Hamiltonian is not necessarily a good approximation. This confirms the picture 
that for small fields, the B-orbits are stable for small F and unstable for small B, with a 
linear transition line between both regimes. 



4.4 The structure of phase space 

The phase space of the averaged Hamiltonian is four-dimensional, and depends on the 
parameters A, F and B. The discussion in the previous sections suggests that the global 
structure of phase space will mainly depend on the ratio %\ we expect that increasing 
both fields, while keeping their ratio constant, will mainly result in an increase of the size 
of chaotic components, without changing the location and stability of the main periodic 
orbits. 

The equations of motion are given by (j§5|), where the expressions of the linear and 
quadratic parts in B are deduced from (fT8|) and ( |88| ) respectively. The manifold of constant 
energy is three-dimensional, and can be represented by a Poincare section. The structure 
of the equations of motion shows that, at least when ^ is not too large, a surface of section 
of the form w e = constant will be a good choice. 



Fig. 10 shows Poincare sections taken at w e = and 



, , in a case with ^ 



1 (sections 

with a smaller ratio of -p look similar). This section should be compared with Fig. |6|a. 
The hyperbolic points located at (wk,K) = (§,0) and (^,0) correspond to the B-orbits. 
Their stable and unstable manifolds separate phase space into two regions, corresponding 
respectively to oscillations around the elliptic BF-orbits and around the S-orbits. In fact, 
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(a) K 



(b) K 




w K 



ir 2jt jt 2jt 

Figure 10. Poincare sections of the singly- averaged Hamiltonian, (a) at w c = and (b) 
at w e = ?, for B = F = 0.1, A = 1 and (if) a — ifo = —0.15 (which corresponds roughly 
to J c ~ 0.5). Hyperbolic points are the B-orbits, elliptic points the BF-orbits. 





FIGURE 11. BF-orbits, obtained by integrating the equations of motion of the unaveraged 
Hamiltonian, (a) f or F = B = 0.05 and (b,c) for F = 0.05 and B = 0.1. Initial conditions 
are given by (106) with e = 0.2. In (a) and (b), F is vertical and B points out of the 
plane, (c) shows a projection on the (B, F)-plane. 



we expect chaotic motions to show up near the separatrices, but they occupy a very small 
area for these values of the fields. 

In contrast to the B-orbits, the BF-orbits move in the (wk, -K")-plane as w e varies. This 
effect did not show up in the doubly-averaged approximation, and is mainly due to terms 
linear in B of the Hamiltonian. Fig. 1C shows that the BF-orbits are roughly located at 
K ~ Qcosz/; e , wk — /SsintUe and K ~ —acosw e , wk — 7r — /^sin-We, where a and f3 are 
of order y (they can be estimated by inserting Fourier series in the equations of motion) . 

In order to understand the geometry of BF-orbits, we observe that 

1. if w c = 0, K reaches its maximum and wk = 0; we know by (Q) that in this case 
the eccentricity is maximal and cos a; = 0; since K ^ 0, however, e is strictly smaller 
than 1, and thus the orbit no longer approaches the nucleus, as it does for B = 0; ( |80| ) 
shows that f2 = ^ or ^ (depending on the sign of J e ); this means that the plane of 
the Kepler ellipse contains B A F, and since cos a; = 0, its major axis is in the plane 
(B,F). 

2. if w e = ^, K = and wk = we also have cos a; = 0, but this time the eccentricity is 
minimal; (]80| ) shows that sinfi = 0(f3), meaning that the plane of the ellipse contains 
F, but is slightly rotated with respect to the (B, F)-plane, by an amount of order -p . 

In fact, when -p is small, we can deduce from (|(]) that O is close to or ir for most values 
of w e . It approaches a step function of the form \ [1 — sign(sin?x; e )] as % tends to zero, 
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Figure 12. Poincare sections at w c — 5 of the averaged system, represented on the sphere 
k\ + «2 + k§ = (A — | J |) 2 for increasing values of y, compare Fig. ||. In all cases, A = 1, 
(H) A = -| - |F, and B = 0.1. (a) F = 0.1, (b) F = 0.05, (c) F = 0.038 and (d) 
F = 0.03. Points are sparse near the poles, as we only plot points obtained with sufficient 
numerical accuracy. 



which means that the orbit approaches the (B,F)-plane in this limit. For increasing B, 
however, the orbit gains some thickness in the direction perpendicular to the plane, and 
rotates around B (Fig. 11). BF-orbits are thus truly non-planar when B > 0. 

Initial conditions producing BF-orbits can be constructed in the following way. Pick 
an eccentricity e G [0, 1). Starting with w c = 0, we have K = and $7 = wk = P can 
be read off the Poincare section (though, the orbit being elliptic, starting with a slightly 
wrong wk will not have dramatic consequences). Taking M = as initial position on the 
ellipse, the initial conditions are then obtained from (|TT|) to be 



x 

y 

z 






A 2 (l 



Px 

Py 
Pz 



COS ft-fc 

sin Qjr 




0. 



(106) 



We can also start with w e = 0, M = 0, and read K = a off the Poincare section, and take 
as initial conditions 



x 







A0Y 2 (l-e 2 )- K 2 J±& 



Px = - 

P y=0 

Pz=0. 



(107) 



The boundaries of the section are given by the condition \K\ = A — \ J e \, and depend 
on wk and w e because J e is no longer constant. They correspond to the location of 
the S-orbit. The equations of motion become singular as \K\ — > A — \ J e \, which makes 
unreliable the numerical computation of orbits approaching the S-orbit. This singularity 
can be tamed by a canonical transformation (J e , K;w e ,WK) *— > (</ e , J ; <fi, —wk)- In the 
case K — > A — J e , J c ^ 0, for instance, it is given by J = A — J e — K and = w e — wk, 
the other cases being similar. We will not elaborate on this point here. 



Fig. |12| and Fig. |13| show phase portraits for increasing values of -p. They mainly differ 
by the behaviour of the B-orbit and its stable and unstable manifolds. While for small 



%, these manifolds connect the B-orbits (Fig. ITOland Fig. 12 a), they connect each B-orbit 
to an S-orbit for larger y (Fig. [E^b and Fig. |13|a,b). When this ratio increases further, 
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Figure 13. Poincare sections of the averaged Hamiltonian for increasing values of 
(a,c,e) are taken at w c = and (b,d,f) at w e = \. Values of A, H and B are the same as 
in Fig. H and (a,b) F = 0.05, (c,d) F = 0.038, (e,f) F = 0.03. The equilibrium points 
K = 0, 10k = ^ and ^ are B-orbits, and the elliptic points near wk = 0,tt are BF-orbits. 
The SB-orbits can be recognized on figures (b), (c) and (d). 
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each manifold folds back to the B-orbit it emerges from, enclosing an elliptic orbit that 
we called SB (Fig. |i~2] c and Fig. [i~3|c,d). Finally, the SB-orbits bifurcate with the B-orbit, 



which becomes elliptic, as we already know (Fig. 12 i and Fig. |l3|e,f). This scenario is very 
close to the scenario obtained for the doubly-averaged system, compare Fig. ^. 

The BF-orbit exists and remains elliptic for fairly large values of -p. It is possible that 
it exists for all values of the field, and transforms into the Z-orbit or C-orbit of the Zeeman 
effect as F — > 0. Verifying this conjecture would, however, require an understanding of the 
transition between the Poincare sections at w c = constant, valid for small and moderate 
p, and those at £1' = constant, valid at large p. At least, it is true that the geometry of 
BF-orbits for w c = is compatible with that of the Z-orbits for D,' = 5 . 

The structure of phase space is thus essentially organized around the B-orbits, BF- 
orbits and S-(or SB-)orbits, which are periodic orbits of the averaged Hamiltonian. The 
other orbits of the averaged system are either quasiperiodic with two frequencies, or belong 
to resonances or chaotic components, which are small, however, for the parameter values 
compatible with the averaging approximation. Intermediate values of p are the most 
favorable for diffusion in phase space, as the heteroclinic connections between B- and 
S-orbits allow fast transitions between various regions of phase space. 

The B-orbits are called S + and S- by []FlWe96 ] , who also describe a periodic orbit 
located in the plane perpendicular to F, called S±. We did not find such an orbit (which 
would correspond to J e = and w c = or ir) in the parameter range we investigated. It 
might be that the BF-orbit behaves as the Sj_-orbit for F <C B. 

Note that going from the averaged to the unaveraged dynamics will add a time scale. 
Hence, periodic orbits of the averaged system may correspond to quasiperiodic orbits of 
the unaveraged Hamiltonian with two frequencies, or to soft chaotic components. Some of 
the quasiperiodic orbits of the averaged system will support KAM-tori of the full system, 
which correspond to a quasiperiodic motion with three frequencies. 



5 Conclusions and outlook 

Action-angle variables and the technique of averaging help to understand the dynamics 
of the hydrogen atom in crossed electric and magnetic fields in two ways. 

The first way is in terms of adiabatic invariants. The semi- major axis A 2 of the Kepler 
ellipse is a constant of motion of the averaged Hamiltonian, and thus an adiabatic invariant 
of the full Hamiltonian on the time scale B~ 2 , compare (|S2"|). The same role is played by 
the averaged perturbing function (H)\ — Hq. In the case B <C F, we have the additional 
invariants J e and ((H)) (see (|S8|)), which evolve on the time scale B^ 1 . In the case 
F -C B, K' = ~L ■'S/B and the averaged quadratic Zeeman term are adiabatic invariants 
on the time scale F^ 1 . 

The second way is in terms of periodic orbits of the averaged system, which organize 
the structure of phase space. The B-orbits, contained in the plane perpendicular to B, 
exist for all fields. They are unstable for B — > and stable for F — > 0, with a roughly 
linear transition line in the (F, -B)-plane. In the limit B — > 0, we also find the S-orbits 
of the Stark effect, and the stable BF-orbits contained in the plane of B and F. As B 
increases, S-orbits appear to become unstable by expelling a pair of "SB-orbits", which 
are then absorbed by the B-orbits. BF-orbits are non-planar for B > 0. They exist and 
are stable in a large domain of field values. It is unclear whether they transform into one 
of the periodic orbits of the pure Zeeman effect as F — > 0. 
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Diffusion in phase space is most prominent in the regime where F and B are of com- 
parable magnitude, when neighbourhoods of the S-orbits are connected by heteroclinic 
orbits. In the other regions of phase space, close to the stable orbits, orbits of the aver- 
aged system are trapped inside KAM-tori. In the real three-degrees-of-freedom system, 
however, diffusion becomes possible in these regions as well, although on a longer time 
scale. This aspect of dynamics still has to be better understood, as do the quantum 



and semiclassical mechanics of the problem which continue to be of interest [vMU97b, 
vMFU97b| , |Nfc9l , |T97| , |CZRT99fl . 
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A Averaging 



A standard result on averaging [V96] is 
Theorem A.l. Consider the initial value problem 

x = ef(x, t) + e 2 g(x, t, e), x(0) = xo, (108) 
where f(x,t) is 2ir -periodic in t, and x,xq 6l". Define the averaged system 



y = e(f)(y), y(0) = x , 



(f)(y) 



2tt 



2tt 



f(y,t)dt. 



(109) 



// / and g are sufficiently smooth and bounded, then 

1. x(t) — y(t) = 0(e) on the time scale 1/e; 

2. Ifyo is a nondegenerate equilibrium of ( 109 ), then ( p.08| ) possesses an isolated periodic 
orbit 7(t) = yo + 0(e), with the same stability as yo if yo is hyperbolic. 



Consider now a Hamiltonian of the form 

H(I 1 ,...,I n ;tpi,...,cp n ) = ilo(-fl) +sHi(I 1 ,... ,I n ]tpx,.. 
where the <pi are angle variables. If we introduce the Poisson bracket 

df dg df dg 



i¥>n), 



-^Ldtpidli dhdipi 

the equations of motion can be written in the form 

ip 1 = H^h) + s{i Pl ;H 1 } 
<Pj = E{(pj;Hx} 
i l = e{I l ;H 1 } 



(110) 



(111) 



j=2, 



,n 
, n. 



(112) 
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If {<pi; Hi} is bounded, H (I\) ^ and e is small enough, we may reparametrize the orbits 
by (pi instead of t, giving 



dL_ = {J i; gi} 

d^ 6 ^(J X ) 



J = 2, . . . , n 
i = 1, ... ,n. 



(113) 



According to the theorem, the dynamics of this system are well approximated by those of 
the system averaged over the fast variable (p\, 



d<fj = Wj] (Hi)} 
M H> (h) 
_ £ {I i ;(H 1 )} 



dtp! 



H (h) 



j = 2, . . . , n 



, n, 



(114) 



where we have introduced 



1 f 2n 

(H 1 )(I 1 ,...,I n ;ip2,...,(p n ) = — / H 1 (I 1 ,...,I n ;ip 1 ,...,ip n )d(p 

^ Jo 



(115) 



We now observe that the canonical equations associated with the averaged Hamiltonian 
(H) = H + e(H 1 ), 



tp l =H' Q {h)+E{y x ;(H x )} 
h = 



j = 2, . . . n 
j = 2, . . . n 



(116) 



are close to 0(e 2 ) to those of the averaged system (114). Thus the dynamics of the 
averaged Hamiltonian are a good approximation, in the sense of the averaging theorem, 
of those of the initial Hamiltonian. In particular, the averaged Hamiltonian (115) is an 
adiabatic invariant of the initial system. 

The averaging procedure can be extended to higher orders in e by the method of 
Lie-Deprit series p69) , p70|| . 



B Poisson brackets 

In this appendix, we discuss the computation of Poisson brackets involving various sets of 
variables used in this paper. 

Since Hamiltonians are usually expressed in terms of auxiliary variables, such as the 
eccentricity e or inclination i, it is useful to start by computing some Poisson brackets 
involving these quantities. Some of them are shown in Table ^. 

For instance, for the perturbation term H\ = j^r 2 [l + cos 2 i + sin 2 i cos(2u; + 2v)] of 
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Table 3. A few useful Poisson brackets involving Delaunay variables. We show Poisson 
brackets between columns and lines, for instance the upper right element is {u; e}. Brackets 
involving G, K and £1 are trivial to compute. 



the Zeeman Hamiltonian, we find 

{A;#i} = 



- A 4 e [l + cos 2 i + sin 2 i cos(2w + 2v)] sin E + -A 3 G sin 2 i sin(2w + 2i 



1 



1 /2r G 2 \ 
{M;Hi} =-r[-r - —cosE) [l + cos 2 i + sin 2 i cos(2w + 2v) 



{G; Hi} =-r 2 sin 2 i sin(2o; + 2v) 
8 

1 /It 
-.—r — 
8 VA Ae 

1 rG 

sin 2 i sin(2w + 2v) sin E 

8 e 

1 rG 

{co; Hi} =— — [l + cos 2 i + sin 2 % cos(2u> + 2v)] cos E 
8 e 

1 r A 1 r 2 

H sin 2 i sin(2w + 2t>) sin£ + -— [-1 + cos(2w + 2v)l . 

8 e 8G 



(117) 



There is an apparent singularity at e = 0. It can be removed, however, by introducing 
variables J = A - G = 0{e 2 ) and (ft = M + to. The transformation (A, G, K; M, u, Q) i-> 
(A, J, K; (ft, —uj, $7) is canonical. One can check that (ft is finite in the limit e — > 0, and that 
the variable £ = Je lw satisfies £ = C( e )- This shows that the circular orbit is indeed a 
periodic orbit of the Zeeman Hamiltonian, as suggested by the averaged system. 

Another way to deal with the singularity of Delaunay variables at e = in the averaged 
case is to use coordinates (^1,^2,^3) introduced by [JCDMW87], see equation (24). Table 
HI gives some useful Poisson brackets involving these variables. 



C Angle variables to first order in F 

For F > 0, the extremal values of £ 2 are given by the condition 

K 2 - 2ai(F)£ 2 - 2H(F)^ + 2F£ 6 = 0, 



(118) 
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Table 4. Some Poisson brackets involving the variables & defined in (| 

where the quantities 

#(F) = -^-3FAJ e + 0(F 2 ) 

«1,2( 



( F ) = 2J ^ +E: ± FA 2 [6 J 5 J, + 3K( J 5 + J„) + K 2 ] + 0(F 2 ) 



(119) 



are obtained by solving equations (43) perturbatively. We find that £ 2 varies between 
limits a\ (F) ± b\ (F) given by 



a x (F) = ai + ±FA 4 [J 2 + 4AJ e - 5A 2 + if 2 ] + 0(F 2 ) 
h(F) =h[l- ±FA 3 (5A + J e )] + 0(F 2 ). 



(120) 



Likewise, the bounded orbits of rf vary between limits (12(F) ± 62(F) which are obtained 
by changing the signs of J e and F in ( |120| ). We may parametrize the level curves of ax t 2 
by 



£ 2 = ai(F)- 61(F) cos V 
r/ 2 = 02(F) - 62(F) cos x, 



(121) 



with the momenta given by (42). The derivatives of the action (54) can then be computed 
as before, with the result 



W A 



: (bi simp + 6 2 sinx) 



i/V 



+ x _ 

"2 2A 2 ' 

+ \F [A(13A + 7 J c )6i sin V - 6 2 sin 2^ - A(13A - 7J C )6 2 sin x + h\ sin 2%] 
+ 0(F 2 ) 

x~4> 



(122) 



+ FA 2 [61 sin if> + 62 sin x] + 0(F 2 
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+ 0(F 2 ). 
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